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Abstract. The aim of this paper is to show that the invariant measure for 
a class of one dimensional chaotic maps, T(x), is an extended solution of the 
Schroder functional equation, q(T(x)) = Xq(x), induced by them. Hence, we give 
an unified treatment of a collection of exactly solved examples worked out in 
the current literature. In particular, we show that these examples belongs to a 
class of functions introduced by Mira, (see text). Moreover, as a new example, we 
compute the invariant densities for a class of rational maps having the Weierstrass 
p functions as an invariant one. Also, we study the relation between that equation 
and the well known Frobenius-Perron and Koopman's operators. 
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1. Introduction 

The notion of an invariant measure plays a central role in the statistical 
characterization of chaotic dynamical systems. In spite of the well established 
theorems about the existence of such measures for abstract dynamical systems the 
construction of such measure for a given dynamical system is in general a difficult 
task. In the case of chaotic maps it is well known that if the invariant measure has 
a density this one is a solution of the Frobenius-Perron functional equation. In one 
dimensional dynamics, some rigorous results about the existence of such measures 
are known, see for example @] EJ HU [TO HOI [HI H2] • There is not doubt about the 
physical significance of the invariant measures, with them we can compute ensambles 
averages, correlations functions, etc., in particular modeling non-equilibrium states [5]. 
Therefore to find reliable analytical methods to construct these measures has special 
interest. 

In the last three decades a series of papers, devoted to deterministic chaos, 
are reporting examples of chaotic maps that are handled by means of a change 
of variables method 0[Tl[Tg[M[]]2[IB[22[2^ Their study 

concerns with nonlinear piecewise transformations that are topologically conjugated, 
or semi-conjugated, to linear piecewise maps. Hence, the invariant density for the 
nonlinear case is obtained by transforming the constant piecewise density of the linear 
one. In particular, it has been remarked in some of these papers that the invariant 
measure of the nonlinear map is in fact the function doing such conjugation. Also, 
this procedure allows to generalize some results proved for piecewise linear maps 
to nonlinear maps conjugated to them. For example, existence theorems, ergodic 
and asymptotic properties or spectral decompositions of nonlinear transformations 

mug Eg. 

Moreover, from the experimental point of view, one-dimensional maps are useful 
to modeling nonlinear processes in diverse fields of science. Also, from theoretical 
point of view with the study of nonlinear maps, which are conjugated to piecewise 
maps, we gain insight to understand integrable chaos p~3] [23] . 

In this paper we provide a new approach to the problem of computing invariant 
densities and measures for chaotic maps, based on Schroder's functional equation. 
This method becomes an alternative to the traditional one given by the Frobenius- 
Perron equation, but also more important, that can be interpreted as an adjoint 
problem to this one. The main difference between these two approaches comes from 
the fact that, for Schroder's equation we only need to know the transformation and 
its derivative, without the explicit knowledge of the inverse transformation and its 
respective derivative, as it is demanded for the Frobenius-Perron equation. Also, it is 
very interesting to us the fact that Schroder's equation is older than the Frobenius- 
Perron equation. Then, we are giving a new look to an old equation. 

This paper is divided in seven parts, section II is concerning with a brief exposition 
of the well known conjugation property between maps. The section III is presenting 
several exactly solvable examples. The fourth one, presenting the main result: the 
formal relation between the Schroder and Frobenius-Perron equations, hence, we 
proceed to the construction of the invariant density. In section V we give a new 
example, a class of functions having the Weierstrass p{x) function as an invariant. 
Finally, we compute the Lyapunov's exponent of these maps and discuss some possible 
generalizations of our method. 
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2. Conjugate maps 

Let us consider a transformation T of an interval / into itself, which preserves the 
measure fi, i.e. ^(A) — u(T~ 1 (A)). Let x e I, the iteration 

x n+1 = T(x n ) , (1) 

defines a discrete dynamical system. Now, if there is a transformation S : I —> I related 
to T by a change of variables T(h(x)) — h(S(x)), where h : I — > / is a continuous 
one-to-one function, we say that the maps T and S 1 are topologically conjugate. If 
in addition, h is a piecewise differentiable function then their invariant densities are 
related by: 

PT{X) = \h>oh-i(x)\ ■ (2) 
The simplest examples of linear piecewise transformations of an interval are: the 
Renyi transformation R r (x) := {rx} and the piecewise continuous transformations 

HB mi [21] 

N r (x) := (-l)^{rx} , (3) 

where [rx] is the integer part of rx and {rx} means rx mod 1 = rx — n, where n is 
the largest integer such that rx — n > for any r S N. 

In the following we are concerning with the case when S is R r or N r , hence 
Ps(x) — l[o,i](a;), for any r e N. Therefore, pr(x) — 1/ \h' o h~ 1 (x)\. Also, we denote 
by T r any map conjugated to N r or R r . It is implied by the conjugation property that 
each map T r has r monotonic pieces that map the interval I onto itself. Such that, 
the inverse function T^ 1 has r monotonic branches on /, each one denoted by TT ■ 

3. Exactly solvable maps 

To show the essence of our method we work out in detail the logistic transformation. 
It is well known that it is conjugated to N2 and that its invariant density satisfies the 
Frobcnius-Pcrron equation. 

3.1. Logistic map 

Let us consider the transformation F{x) = 4x(l — x), where x E [0, 1]. The function 
h{6) = sm 2 (f 9), 9 € [0, 1] satisfies: 

F[h{9)] = h[N 2 {9% (4) 

hence, by flU) 

P p ^ = /\ ( 5 ) 

7rya;(l — x) 

Finally, the invariant measure for F is 

r 1 

^( a; ) = / — / „ . dy, x e I, (6) 
Jo TT^y{l-y) 

2 

= — arcsin(v / x) , (7) 

7T 

in other words, /j,(x) — h^ 1 (x). It should be noted that (HJ can be written as 
F(x) = h[N2(h (x))] then, by periodicity of the sin function the last expression 
is equal to F(x) — h(2h^ 1 (x)), i.e. F(x) = sin 2 (2sin _1 y/x). 
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3. 2. A collection of exactly solvable maps 

Exactly solvable iterations are known since the nineteenth century, the interest was to 
understand Newton's method [SI GO HI]- However, in the last thirty years a boom on 
the search of analytic results for chaotic systems shifted the attention to the topic. 

A very important examples are provided by the Chebyshev polynomials. It is 
well known that they are ergodic and mixing maps [22j . From its definition: T r {x) = 
cos(r cos" 1 x), we can see that they are conjugated to N r (x) maps. To show this, we 
start with the function h{0) = cos(7r6'), and obtaining T r (x) = h[N r (h" 1 (x))]. Hence, 
using ([2]) pt{x) = l/ir^/l — x 2 . Such that, its invariant measure is p(x) = — arccos(a;). 
We can see also that F(x) is also conjugated to T 2 {x) via N 2 (x). 

Now, we list some papers on the subject: a) in |15j the invariant density of 
maps conjugated with the Jacobian elliptic sine function, sn{x) 1 are presented. The 
class of maps studied is T r (x) = f{sn{rsnr 1 {f^ 1 {x)))), where f(x) is an appropriate 
function mapping the interval [0, 1] onto itself; b) in [19) some examples of one 
dimensional iterations are solved using the addition properties of trigonometric and 
elliptic functions; c) in [14] transformations of the real line associated to special 
statistical distributions arc built. In particular, the Cauchy, F and Z distributions are 
associated to maps of the Real Line conjugated to N 2 ; d) in [BJ [28] Newton's algorithm 
to search the square root of \/^3 is found to be conjugated to N 2 , or as in [57] to 
2x + 1/2 mod 1; e) recently, in [29] a similar method as [15] is given; f) finally, an 
exception is the paper [25] where Schroder's equation (see below) is explicitly used 
in order to generalize the properties of the Chebyshev polynomials to more general 
functions. But, in that paper the relation to the Frobenius-Perron equation is only 
suggested. Nevertheless, it should be noted that the invariant densities of all those 
examples are exactly computed exploiting the addition and periodicity properties of 
trigonometric and Jacobian elliptic functions jointly with the conjugation property 
with R r or N r . It is worth to mention, that the examples handled in these papers 
were obtained by each author mostly in an independent way to the others ones. For 
that reason we speak about rediscovered maps in the last paragraph of Section [S] 



4. Schroder's equation 

By definition, if p(x) has a density p{x), then 

K x ) = / p{y)dy , (8) 



Now, the preservation of p under T implies that p(x) satisfies the Frobenius-Perron 
functional equation [llj : 

P (x) = ± , ggQjgl , (9) 

or, in an integral form: 

p(x) = f 6(x-T(y))p(y)dy. (10) 



Now, using this last expression, we ask the question about the left eigenvalue problem 
for the Frobenius-Perron equation, which is defined by: 

A<?(z) = / q(y)S(y - T(x))dy = q(T(x)) . 
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This equation is nothing but Schroder's functional equation [51 [T5] : 

Xq(x) = q(T(x)) , (11) 

where A € R and q(x) is a function of a real variable. Hence, the iteration fTJ) is 
linearized, such that its n-th iterate is (at least locally): X n q(x) = q(T n (x)). 

The classical result in the local analysis of solutions of Schroder's equation is 
Koenig's theorem which guarantees the existence of a solution in the neighborhood of 
a fixed stable point of T [XT] US]- However, more difficult is our present interest in the 
global study of that equation, see [U [T71 [18] for a discussion. 

Now we will consider a differentiable map T : I — > I conjugated to N r (x), then 
T has r monotonic pieces on /, and T -1 has r branches on /. We proceed as follows: 
taking the derivative on both sides of Schroder's equation (fTTjl and after taking its 
absolute value, we obtain 

\\\\q'(x)\ = \T>(x)\\q'(T(x))\ . (12) 

Indeed, substituting the j-th branch of in ([T2"|) 

\X\\q , (Tr\ x ))\ = \T , {Tf\ x ))\\q'( x )\ . (13) 

x E I. We observe that this expression is valid for any j = 1, . . . , r. We are now ready 
to introduce the result of this paper, putting together all terms of the form ([TU]) : 



|A| 



q'(T-\x)) 



T'{T-\x)) 



r\q'(x)\, (14) 



where j runs over all the inverse branches of T. As we can see, this is the Frobenius- 
Perron equation for |g'(x)|, therefore, if there exist a function a(x) = \q'(x)\ as a global 
solution of (|14|) . defined on /, for |A| = r, then the corresponding invariant measure 
of T is 

fi(x)= [ \q'(y)\dy, yel. (15) 







It is very easy to prove, by direct substitution, that the examples handled in 
P [III HS1 EE1 [SI H21 123 HZl [2H1 , which are conjugated to N r , satisfies the equation 
(fl~2j) . For example, we use the logistic map to show the method. First, substituting 
F^{x) = (1 ± \/r^)/2 in the lhs of ([12]), 

\Mp{FZ\x)) = |A ' = , (16) 

^F ± \x){l-F ± \x)) ^ 

and in rhs of (fT2l). 



\F'(F£\x))\ p{x) = 4 ^!__ = , (17) 

such that, we have an equality if and only if |A| = 2. 

It is worth mention that the ([T2]) is formally different from ITT]) . However, if the 
functional equation a(T(x)) = jT(x) a i x ) nas a continuous solution a(x) in /, fulfilling 
A = T'(x), then q(x) = a(x)dx is a C 1 solution of pT]) . Thus, there exists a 
one-to-one correspondence between the C 1 solutions of Schroder's equation and the 
continuous solutions of (fT2l) 1181. 
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5. Rational transformations on the Real Line 



Now, as a non trivial example we consider a class of rational transformations having 
the Weierstrass p(z) elliptic function as an invariant one. This function is defined by 
an infinite series, also it is known its duplication formula. 
The Weierstrass p(z) elliptic function is defined as pQ: 

1 



(z — 2mu>i — 2ncL>2) — (2mu>i + 2nuj2) ' 

where z € C and u>±, 0J2 are two numbers the ratio of which is not real, and the 
summation is take over all m,n £ Z but excepting the case when m = n = 
simultaneously. This is an example of a double periodic function with fundamental 
periods 2wi and 2^2- As it is standard, that function is parametrized with two 
numbers, g 2 and 53 which are called the elliptic invariants (and are functions of u>i, 
We describe Schroder's method to obtain a rational function invariant under p{x): 
there is a duplication formula [T]: 



p(2z) = -2p{z) + 



2p'(z) 



such that, using the identities: p"{z) = 6p 2 (z) 
33, we can built a two parametric rational function 



f , and 



y) 2 (z) = 4 P 3 (z)- ff2 pW- 



</2 .</:) 



+ f^ + 2 53 z + (f) 2 



4z 3 - g 2 z- 33 
having p(z) as an invariant function, i.e. p{2z) 



(18) 



Tg 2 .g 3 (p(z)). Now, let us consider 



the case 32,53 S K, then p takes real values on the real line. The inverse function for 
p is also known [3]: 



92s - 33 



Now, if the discriminant A = ^(.9! — 53) is positive, then 4s 3 — g 2 s — — has 
three real and distinct roots and p has the period u\ £ R. Also, the above relations 
are unaltered by the choice z — x. Now, to use the Frobcnius-Pcrron equation 
we must compute the inverse function of the T given by (|18p . and its derivative. 
Instead of this, we proceed to compute the solution of the corresponding Schroder's 
equation associated to such T. We can see that in analogy with section 13.11 that 
— Tg 2 ,gi,{ x )i therefore the invariant density is given by 

' (19) 



p(x) 



d -1/ \ 
ax 



V4a 



52a; - 33 



i.e., the eigenfunction of pip induced by Tg 2 ,g 3 and having eigenvalue |A| = 2. 

For example, if 92 = 4 and 33 = 0, we have the rational transformation 
T{x) = (x 2 + l) 2 /4x(x 2 - 1), where x ^ 0, ±1, then T(p(x)) = p(2x). This example 
was given by Lattes in 1918 [5]. Another important example, is given by 

= Ax(l-x){l- k 2 x) 

(1 — K 2 X 2 ) 2 

which satisfies the semiconjugacy relation T(sn 2 (x)) — sn 2 (2x), where sn(x) is the 
Jacobi elliptic sine function. As it is remarked by Milnor [24] . this example given also 
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by Lattes, was in fact studied by Schroder in 1871 [55]. However, we found that it 
was also rediscovered in [TH [JJ3 12H| ■ A different example, studied by Boole in 1872, 
is the linear fractional transformation 

T(x) — ad — be = 1 , c =/= , 

ex + d 

whose Schroder's equation is solved using an iterative method, called by him as 
Laplace's method, which it is exposed in [7]. For this example, our method gives 
p(x) = 1/(1 + x 2 ) as its invariant density, provided \a + d\ > 2. Again, this is a case 
of a rediscovered example [33] . 

6. Lyapunov's exponents 

Finally, the Lyapunov exponent A(T) of the maps T r is computed using Schroder's 
equation induced by it. By definition, A(T) :— J T In |T/(x)| p(x)dx, then taking 
logarithm of (|12|) . and using the invariance of /i under T r , we have from the previous 
equation the expected result A(T) = lnr. In other words, A(T) is equal to the 
logarithm of the number of monotonic pieces of T r (x). 

7. Discussion and Conclusions 

To the best of our knowledge, we are the first to point out the significance of 
Schroder's equation to compute invariant densities and measures for chaotic maps. 
Our approach makes the main difference with the corresponding work of the authors 
inPOjaiigQIlIigililC^ As it is mentioned in Section[3J2 

the work by Mira [25] using Schroder's equation, is focused in the generalization 
of the properties of Chebysev maps. It should be noted the interest of authors 
working on the subject of functional equations on the problem of computing invariant 
densities, see for example [4]ll7 [ l2T j . On the other hand, authors working in dynamical 
chaos are interested in methods to solve the Frobenius Perron functional equation 
El HH H31 HH HSl [HI [HI EUl EH EH ESI EH EH EH El]. However, in the current 
literature the problem of computing invariant densities or measures appears unrelated 
to the problem of solving Schroder's functional equation. 

We remark that the measure fi(x), which is the integral of |<z'(x)| in ([T5]) , could 
be not equal to q(x), is in that sense that we speak about an extended solution of 
Schroder's equation. 

Also, we note that in Ergodic Theory Koopman's operator is introduced as 
the formal adjoint of the Frobenius-Perron operator on densities in the L°° sense 
[UJ [22l [23] . I n our context, we are considering Schroder's equation as a left eigenvalue 
problem associated to the Frobenius-Perron functional equation, such that our case is 
a more general one. It should be pointed out that this investigation is not limited to 
the one dimensional case. However, higher dimensional dynamics seems to be more 
difficult and may be requiring a different approach, mainly due to the formal definition 
of Schroder's equation in several variables [P2] . 

In this paper we are mainly concerned with real transformations on the interval, 
but another interesting aspect of the subject is the study of the complex Schroder's 
equation associated to chaotic transformations of the complex plane [12l [27] . 

Finally, in the one dimensional case becomes very interesting the study of strange 
attractors at 2™ banded chaos of the logistic map F a (x) = ax(l — x). Here a S (aoo,4), 
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and a x is the accumulation point of the period doubling bifurcation [TBI HH HB]. 
Further studies in this direction will be published elsewhere. 
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